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1 Introduction 



The seminal and, so far, most used way to formulate the quantum theory of a particle 
or a field makes wide use of the Hamiltonian description of classical mechanics |]l|]. The 
standard rules for constructing the momenta and the Hamiltonian function, however, 
cannot be applied when the Lagrangian is singular. In such a case it is not possible 
to extract the functional dependence of all the velocities on the momenta in order 
to obtain a Hamiltonian function of coordinates and momenta only. Dirac's method 
concerns the study of classical systems using the Hamiltonian method when the usual 
procedure fails due to the singularity of the Lagrangian 0]. Dirac gave very general 
rules to construct the Hamiltonian and calculate sensible brackets that can be used 
to describe the classical and, by the canonical quantization procedure, the quantum 
dynamics. 

One of the most interesting situation where Dirac's method of handling singular 
Lagrangians can be applied is in confining particles on curved manifolds P, Part 
of this interest is certainly due to the presence of additional terms which arise in 
many quantization procedures on curved manifolds [0, 01 and is far from being 
clarified. In this letter we will focus our attention on the connection between the 
additional terms which occur in the quantum Hamiltonian and the problem of the 
operator ordering prescription. 

In Section 2 we briefly review Dirac's method of handling singular Lagrangians. 
In Section 3 we quantize a free particle constrained on a circle following the standard 
procedure, i.e. reducing from the very beginning the number of degrees of freedom. 
Then we solve the same (classical) problem using Dirac's method, recovering a new 
set of canonical brackets. Finally we quantize using this bracket algebra, by focusing 
our attention on the construction of coordinates, linear momenta, angular momentum 
and Hamiltonian operators and on related ordering problems and we will finally write 
the Schrodinger equation. Section 4 contains our conclusions. 

2 The Dirac method 

Let us start by outlining the Dirac method [0 and introduce notation. Take a consis- 
tent Lagrangian L{x,x) with coordinates. The classical dynamics is obtained by 
the least action principle: 




■tl 



S[x] 



5S[x\ = 0, 




which in terms of the Lagrangian gives Euler-Lagrange equations 



dL_ 

dxi 



d dL 

dt dxi 



0. 
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We define momenta and Hamiltonian and obtain the usual (Poisson) brackets between 
momenta and coordinates: 



dL 

P^ = g-, {^ = l,...,N) (2.3) 

H{x,p) = '^piXi{x,p) - L{x,x{x,p)), (2.4) 

i 

[xi,pj] = Sij, (2.5) 

and for any function A of x's and p's (not explicitly dependent on time), 

A = [AH]. (2.6) 

Two scenarios are possible. In the typical case one can invert pi{x,x) to obtain 
Xi{x,p)] if this is not possible, not even locally, the Lagrangian is said singular and 
its Hessian with respect to the velocities vanishes 



dL 



0. (2.7) 



dxidxj 

In such a case we act differently. We consider those relations in ( |2.3|) which hinder 
the inversion (this step will be clarified in the example of Section 3) as a series of 
constraints 

^ (2.8) 

which must be satisfied "weakly" (namely, their Poisson bracket with any given quan- 
tity may not vanish) along the physical trajectory. In this way we obtain a number 
(say M) of constraints which Dirac called primary because of their direct derivation 
from the Lagrangian. Notice that a Hamiltonian is required to be independent of the 
velocities. If we are not able to erase the x dependence, then the straightforward ap- 
plication of the hamiltonian method is impossible. To solve our problem we proceed 
as follows. We add to H all our primary constraints multiplied by arbitrary functions 
of time Uj, to obtain the total Hamiltonian Hj' 

M 

Ht = H + Y.uj(I)j{x,p). (2.9) 
i=i 

This could seem to imply an arbitrariness (additional freedoms are introduced) but 
we require a number of consistency conditions: each constraint must be zero during 
the whole evolution, if it is initially zero: 



Ht]^0 (j = 1,...,M). (2.10) 



If these equations are consistent, three cases are possible: an equation can give an 
identity; it can give a linear equation for the uj; it can give an equation containing only 
p's and x's, in which case it must be considered as another constraint. The constraints 
that arise from this procedure will be called secondary, for obvious reasons. Even for 
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these, we impose consistency conditions and this procedure is continued until we have 
a set of identities and hnear equations for the m's. Now we have enlarged our set of 
constraints to include the secondary ones and we have a new number of constraints, 
say K. 

We have by now defined a constraint as a quantity which satisfies 

0, ^ 0, (2.11) 
[(pj^Hr] ^ 0. (2.12) 

This defines a linear vector space (due to the linearity of the Poisson brackets) and 
so any linear combination of constraints is again a constraint. It is of great impor- 
tance for our purposes the distinction between first class and second class constraints. 
The first are defined as the constraints which "commute" (i.e. have vanishing Pois- 
son brackets) with all the other constraints. The second ones have at least one non 
vanishing bracket with some other constraint. It may happen that we can take lin- 
ear combinations of second class constraints and obtain some first class constraints. 
This situation brings to light the presence of some gauge degrees of freedom. Dirac 
showed the profound difference between this two classes. In fact we can switch to new 
canonical brackets in order to set all of our second class constraints strongly equal 
to zero. This means that in any given quantity, such as the Hamiltonian, we can set 
them to zero "by hand". The first class ones, however, will "survive" (even in the 
Hamiltonian with their arbitrary multiplicative functions u). In the following analysis 
we will not deal with first class constraints and so will not discuss them any further. 
Every constraint that we will find will be of the second class. In such a case, we can 
safely change to the new canonical brackets, the so called Dirac brackets, defined as 
follows: let 

M,, = [0„0,] (2.13) 

and its inverse 

G., = {M-% (2.14) 

(the invertibility of M is a particular feature of the absence of first class constraints: 
in general M is defined on the subspace of second class constraints only). Then for 
any two quantities A and B we define the Dirac bracket: 

[A, B]d = [A, 5] - E [A, <P,]G,A<Pi, B]. (2.15) 

These brackets have all the properties of the Poisson bracket plus one: for any dy- 
namical variable A we have 

[A0.]d = O, (2.16) 

A=[A,Ht\^[A,Ht\d, (2.17) 



as is easy to see (for ( p.l7| ) use ( p.l2| )) 



The very meaning of this redefinition of the canonical brackets is simply a change 
of variables from the original phase space to the constrained manyfold H. Having 
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obtained a set of canonical brackets, we can now quantize, by looking for self-adjoint 
operators which satisfy the canonical commutation relation (each quantity in the 
righthand side must be multiplied by ih). 
Let us now look at an interesting example. 



3 Particle on a circle 

3.1 The standard approach 

We want to quantize the following free particle Lagrangian 

L = l{x^ + f), (3.1) 

subject to the relation 

r^ = x^ + y^ = rl (3.2) 

{tq being a positive real constant) which must be satisfied at any time. This describes 
the motion of a particle of unitary mass in the xy-plane, constrained on a circle of 
radius tq. We can make a change of variables, from cartesian to polar coordinates 

X = rcos9, (3.3) 
y = r sin 6, 

after which, using (|3.2|) , the Lagrangian reads 

L = Kle\ (3.4) 

We have now a new Lagrangian with only one degree of freedom 6. We can define 
the momentum pe 

Fit 

Pe = ^ = rle (3.5) 

and the Hamiltonian 

2 

H{e,pe)=epe-L = ^. (3.6) 

The radial degree of freedom r disappears (as implicitly did any other non-dynamical 
degree of freedom, such as the z coordinate in (pTTI)). The Poisson bracket is 

[e,Pe\ = l. (3.7) 

Now, let us quantize: define two self-adjoint operators 9 and pe satisfying the 
canonical commutation relation (CCR) (^=1): 

[e,pe\=t (3.8) 
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(we shall use the same notation for Poisson brackets and commutator of operators, 
since no confusion can arise). We can find such a couple of self-adjoint operators in 
the Hilbert space Ti = L'^{0, 2tt) and their expression is: 

§■^{6) = e-ipie), (3.9) 

We add the constant a in the momentum pg to mimic the possible presence of a 
magnetic field enclosed in the circle (see the discussion after ( |3.44| )). Their domains 
are chosen to be respectively Dg = H and Dp^ = {■?/'€ Hlipi^O) = 4'{27i),ip' G H}. 
These are dense subsets of Ti. Notice also that we have chosen one of the infinite 
self-adjoint extensions of the momentum pg. The Hamiltonian reads 

and is self-adjoint in the domain of pg, i.e. Dp^. The Schrodinger equation is (rein- 
serting m and h) 

This is what we expected. 



3.2 Dirac's approach 

Let analyze the same problem with Dirac's method. We start from classical dynamics. 
We want to find the extremum of the action with the Lagrangian defined in ( p.l|) , 
subject to the constraint 

(j) = x"^ + y"^ - rl ^ 0. (3.12) 

We use the method of Lagrange multipliers and search for the extremum of the 
action with the new Lagrangian 

L{x, X, y, y, A) = + - \ix^ + - r^), (3.13) 

the quantity A being treated as an additional dynamical variable. This Lagrangian 
gives rise to an action functional S[x,y, A] which must be varied with respect to x,y 
and also the "new" degree of freedom A. If we want to use the Hamiltonian method 
with this Lagrangian, we must start by calculating the momenta: 





dL 




Px = 


dx 
dL 


= X 


Py = 


dy 


= V: 




dL 


= 0. 


Px = 




~d\ 





(3.14) 
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It is apparent that we are facing the situation discussed in the Introduction and in 
Sec. 2: one of the momenta disappears. So we proceed as previously sketched: read 
the relation pa ~ as a primary constraint: 

0i=Pa~O. (3.15) 

This is our only primary constraint. Build up the Hamiltonian 

H = p,x + pyy + px\-L = ^ + ^+Px\ + \{x' + y^- r^). (3.16) 
We now include 0i multiplied by an arbitrary function of the time ui: 

1 2 

i^T = ^ + ^ + A(a;2 + - rl) + u^x- (3.17) 

Notice that A has been absorbed in the arbitrary function u\. The consistency con- 
dition ( plO| ) is 

^ 01 = [01, Et\ = [pA, Ht] = -{x' + y'- rl) , (3.18) 

which is a new constraint, that the Lagrange multipliers had already implicitly im- 
posed (0 in ( CT) ) 

02 = </ 

The consistency conditions ( p.lOp for 02 yields 

03 = xp^ + ypy^0 (3.20) 
and by imposing ( |2.10|) also for 03 we obtain 

<p4 = pl+pl-2{x^ + y^)X^0. (3.21) 



x^ + y^ -rl^O. (3.19) 



ui = Y—^{xp, + yPy) ~ 0. (3.22) 



These are additional constraints. If we impose (p.lO|) for 04 we get an equation for 

Ui. 

2A 

x'^ + 1/^ 

Since in the following we shall use only Dirac brackets we regard any constraint as 
a strong equation and drop the term Ui0i from the total Hamiltonian. We can also 
drop the term containing the Lagrangian multiplier because of 02. So our Hamiltonian 
becomes the free one: 

^T = f + f. (3.23) 

The fact that the Hamiltonian function of the constrained dynamics is exactly that 
of an unconstrained dynamics may seem strange. One could (erroneously) argue 
that even the equations of motion would be the same. This is not correct because 
we will change the canonical brackets. All additional information characterising the 
constrained dynamics is now contained in these new canonical brackets. One could say 
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that Dirac's method "drains" information from the Lagrangian, where it is contained 
in the additional degree of freedom A, giving it to the canonical brackets, where it 
is contained in a non-trivial algebra. In this process, however, the information on 
the topology of the problem is made explicit, as we shall see in the short discussion 
just after the algebra construction. This point of view is very useful in quantum 
mechanics. 

We have four constraints and what we need now is the algebra of the Dirac's 
brackets. We calculate the matrix (r = {x,y) and p = {px,Py)) 



M 



and invert it to get 



G 







-2r2 










-2r2 
-4p ■ r 





2r2 


4Ar2 



2p2 



(p2 + 2Ar2)/2r4 
—r ■ p/r^ 
l/2r2 



-(p2 + 2Ar2)/2r4 


l/2r2 




r ■ p/r^ 
-l/2r^ 





2r2 
Ap ■ r 

h4Ar2 




-l/2r2 






(3.24) 



(3.25) 



We can now calculate the Dirac brackets of any two quantities and appreciate their 
physical meaning. 

To start off, let us first consider an interesting example of the difference between 
Poisson and Dirac brackets. We can check whether ( |2.16| ) is true for (pi = px and 
A = X. The commutation rule between the Lagrange multiplier and its momentum 
changes from [X,Px] = 1 to 



[A, Pa] 



D 



1 - Hf^^^iJCijl^J^PA] 



l-[X,4>i]Gu[(p4,Px] = 1-1 



1 
2^ 



-2r' 



0, 



which enables one to see how the Dirac brackets work in order to satisfy the constraints 
strongly. We also find (we have replaced r with ro in each quantity by using 02 = 0): 











1 


[y,Py]D = 


1 


y' 






xy_ 

rl ' 


[y^Px]D = 




xy_ 

rl ' 


[x^y]D = 





^{xpy - ypx 


[Px,Py]D = 





(3.26) 
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This brackets have a nice geometric interpretation. According to the Poisson bracket 
= 1; Px is the generator of translations along the x axis. However this property 
cannot be preserved in the constrained algebra, because typically we cannot translate 
in the x direction while remaining on the circle. This can be done only at the points 
(x = 0, 1/ = ±ro) where the first and fourth equations of ( p.26| ) reduce to the Poisson 
algebra. Another feature is to be noticed: x and y still commute. We can understand 
this because x and y are the generators of translations in the corresponding p's di- 
rections; however there is no constraint containing only momenta so any given point 
in the p^jPy-plane is allowed, by suitably adjusting the other coordinates x^y and A. 
This is not the case of the coordinates x and y, as one can readily see: for example, 
the point x = 2rQ,y = tq is not allowed even by making additional translations of 
momenta and A, because of 02- 

We can write the Hamiltonian in the form ( p.lO| ) defining L^: 



Lz = xpy - yp^. (3.27) 
Squaring it and using 02 we obtain 

= rlipl + pD - {xp, + ypy)', (3.28) 

and using 03 we obtain 

^=^(Px.+P?) = ||. (3.29) 

One can identify with po by writing coordinates and momenta as functions of 6 
and Lz- This can be done by solving the equations ( |3.27] ) and ( |3.20D for p^ and Py 
and using (|3.2| ). We get: 

X = rocos6', y = rQ sin 9, (3.30) 

1 1 

Px = L^siiaO, Py = — L^cosd. (3.31) 

To ro 

The reader can verify that all the relations obtained by the Dirac brackets algebra 
are equivalent to the single bracket [9,Lz] = 1 (e.g. [x,Lz]d = —y should be read 
[cos 9, Lz] = — sin 6^ and so on). 

Equations (p.26|) pave the way to quantization. We shall see that the quantization 



of the Dirac algebra is not a trivial problem: our recipe will be the requirement that 
some operators be self-adjoint (or at least Hermitian). This requirement will play a 
fundamental role in our analysis. We look at an explicit representation of the self- 
adjoint operators x, y,Px,Py (notice that we will not deal with the operators p\ and A 
because they are completely defined by 0i = and 04 = respectively) satisfying this 
algebra. We must, however, impose the (now) strong equalities 0j = (^ = 1, 2, 3, 4). 
So (in the following we will drop all hats on operators), = + = Tq and there 
exists a self-adjoint operator 6 on the Hilbert space Ti. = L'^{0, 2tt) such that: 

X = tqcosO, (3.32) 
y = tq sin 9. 



8 



We will determine the momentum operators in order to satisfy the following equa- 
tions: 



[y,Px] 

[Px.Py] 



-i — 

.xy 



(3.33) 



- — {xpy - yp^ 



Using the fact that {de stands for the ^-derivative, F for any n-times differentiable 
function) [dg, F{6)] contains derivatives of order less than or equal to n — 1 and looking 
at the first two equations in (|3.33|) (whose right hand side does not contain momenta) 
one can infer that the p operators in the 6 representation contain only first order 
derivatives. Then, in the most general case. 



ro ou To 

py = --9{0)^ + -b{e). 

ro oO ro 



(3.34) 
(3.35) 



Using these expressions we solve for the unknown functions /, g, a and b. The first 
equation in (|3.33| ) yields 



d 



d 



cose -tfio)— - -if{,e)— cose = ism'e 



de 



which is solved to give 



de 



sin e. 



Analogously, the solution of the second equation in ( |3.33D gives 

g{e) = cose. 



(3.36) 
(3.37) 
(3.38) 



At this stage the third, fourth and fifth equations in (|3.33|) are identities and yield no 
information on a and b. However, some insight on their form can be obtained from 
the last of ( p.33| ), which gives 



a cos e + b' sin e = —b cos e + a sin e, 

where the primes denotes derivatives. This yields 

a' = -b, 
b' = a. 



(3.39) 
(3.40) 
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However, there are other equations which must be satisfied: 



[x,H] = ip,, (3.41) 
[y,H] = ipy. 

These are hnearly dependent and both equivalent to 

ia cos 6 + ib sine = (3.42) 

By using (|3.41|) this turns into an equation for a whose solutions, under the additional 
requirement that px and Py be hermitian operators (this is a necessary step in order 
to require their self-adjointness), are: 

a{6) = - cos6' + a sin ^, 
b{6) = —a' = -sm6 — acos6, 
where a is an arbitrary real number. Putting all the results together we obtain 

Px = — sm9— + — cos9 -\ sin6', (3.43) 

t d 1 ex 

Py = cos6'— H sin 6* cos^. 

ro oO 2ro 

We can put these equations in a compact form by using the anticommutator (for any 
operators A and B: {A, B} = AB + BA): 

Px = ^e^°^|z^,sin^|e-*"^ (3.44) 

Py ^ ^e-{4,-cos.)e--. 

Written in this form, these equations readily show some properties of these operators. 
First, they are the Weyl ordered operators of the classical quantities (|3.31|) but this 



ordering arises naturally by taking suitable solutions of the algebra equations. Second, 
these p's are self-adjoint in the domain Dp^ defined after ( |3.9|) . Finally, equations 
( |3.44 ) also show that different p's, corresponding to different a's, are connected to 



each other by means of gauge transformations; this interesting property can be easily 
related to the Aharonov-Bohm effect (see |^), identifying a with where $b is 

the flux of the magnetic fleld enclosed in the circle. 

One can check that all the constraints are satisfled: remember that we have chosen 
the expressions of x and y to satisfy 02, set p\ = Oto satisfy 0i, deflned A to satisfy 04, 
so we must manage only with 03. Physically 03 /rg is the radial part of the momentum 
Pr = [r ■ p)/rQ (the vector r being on the circle: = Tq) so we choose to represent 
it with a Hermitian operator. We therefore order it {W stands for 'Weyl ordering' 
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which coincides with any other sufficiently symmetric operator ordering procedure for 
this simple quantity) in order to get the Hermitian expression 



1 



{{x,Px} + {y,Py}) 



(3.45) 



One can easily see, using the solutions 



and (|3.43|), that (^^^w = 0. Con- 



versely, using the algebra relations (|3.33| ) it is possible to show that if a non Weyl- 
ordered expression for 03 is constrained to zero the momentum operators are not 
Hermitian. In fact using the algebra of commutators ( p.33| ) one readily gets three 
equivalent expressions for the momenta operators: 



Px 



^ {-y, L J + —x[-- + xp, + ypy 



1 r r ^ 1 



+ p.j,x + pyy 



Py 



2r§ 

— 2 {-1/,LJ + -2X03,W, 
'^'0 '0 

^ {x, L,} + ;|?/ (-^ + + ypy 

^ {x, L J + Q + p,x + pyy^ 



{x,L,} + ^y4>3,w 



(3.46) 



(3.47) 

'0 

These expression are not Hermitian if we set xp^ + yPy = or p^x + Pyy = 0. On the 
contrary they are Hermitian if we set 03, vi/ = 0. Therefore the hermiticity of these 
operators and their Weyl ordering are strictly correlated. We shall come back to this 
remarkable point in the following. 

We can now build up any quantity we need in our quantum theory, for example 
the z component of the angular momentum 



XPy 



d 



a 



and the Hamiltonian, from (|3.23|) 
1 



H 



-XPx+Pv 



d 



d 



(3.48) 



(3.49) 

\ 2' 



i sin 9-— H — cos 6 + a sin 6 \ + — z cos 6-— H — sin 6 — a cos 



89 



89 



+ 



One can check that the ground state energy is 



2re 



1 

2"" 



+ ^0, 



(3.50) 
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where (in ordinary units) Eq = and a G [0, 1[ , a = a mod 1. Notice that 
we have obtained a constant Eq which was absent both in the classical Hamiltonians 
(|3.29|) and ( |3.6|) and in the quantum Hamiltonian ( |3.1CI|) . Let us now discuss the 



domains in which these operators are self-adjoint. We see from that a good domain 
for its definition is the previously defined Dp^ = {■?/'€ L^(0, 27r)|'0(27r) = ip{0),ip' G 
L^(0, 27r)}. In different domains will not be self-adjoint anymore ad so will not be 
an observable. The Schrodinger equation (reinserting the mass and h) reads 

dtp ( d \^ 

i/i— = 7- 2 ""^^ ~ '-^ i^ + Eoi/j, (3.51) 



dt 2mrn 



which differs from ( |3.11D for the presence of Eq 



At this point one should focus on the connection between the additional term 
Eq and the quantization procedure we have established. First of all, observe that 
if we had not required that the p's be hermitian we would not have found such a 
term. For example, taking H = L^/2rQ on the circle in the classical context and 
then substituting with the corresponding quantum operator gives (of course) no 
additional terms in the Hamiltonian. Another way of getting rid of the additional 
terms is by abelian conversion of the algebra 0, |TT|, by introducing an additional 
set of coordinates. In practice, there are many ways of dropping or changing the 
additional terms without changing the algebra. 

Similar terms arise in the quantization on curved manifolds as an effect of the 
(intrinsic) curvature of the manifold itself, as shown by DeWitt and successively 
elaborated by Schulman |^ . We stress however that the constant Eq found in Dirac's 
procedure cannot be put in direct correspondence with such curvature terms. Indeed 
for the quantization on the circle the scalar curvature is 0, while the additional term in 
Dirac's procedure is I/Stq. However, the fact that no direct proportionality is present 
between these two energies does not mean that one of them (namely Dirac's one) is 
unphysical, but rather that the procedures of quantization leading to them reflects 
different physical processes. While, for example, in the path-integral procedure no 
mention of the embedding space (for example the circle as a subset of the plane) is 
made, in Dirac's procedure this embedding is unavoidable. Moreover a parallel with 
the classical version of Dirac's procedure hints that the physical evolution it suggests 
is like making a small step (free evolution) in the embedding space followed by a 
'projection' on the constraining manifold, obtained by dropping the component of 
the step orthogonal to the manifold p[ . In order to visualize this way of constraining 
the particle on the manifold one should follow the procedure we outlined, i.e. should 
quantize the generators of translations px and Py of the embedding space and then 
calculate the Hamiltonian. On the contrary DeWitt's procedures relies upon the 
presence of a natural metric on the manifold and the additional term is due to the 
intrinsic curvature of the manifold itself. This scenario obviously sets aside a possible 
embedding of our manifold in a bigger one. The choice of a given procedure should 
depend on the physical process one has in mind. 

Additional work is needed in this direction, for the problem is certainly far from 
being solved, as the physical significance of the additional energies appearing in the 
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different procedures is not completely understood. 



4 Conclusions 

As we have shown, the Dirac method yields deep insight even in a simple example like 
the one we considered. The construction of the Dirac algebra of brackets is non-trivial 
and instructive and even more interesting is the search for an explicit representation 
of the self-adjoint operators satisfying the algebra and the constraints. One must look 
at their functional form and identify and interpret any possible freedom inherent to 
their choice. Then one must look at their domains of definition, facing sometimes 
ordering problems. Eventually, one gains a better comprehension of the Hamiltonian 
formalism, the connection between Dirac algebra and the topology of the constrained 
manifold and the quantization procedure on this manifold. An interesting explicit 
result we have obtained is the presence of an additional energy term different from 
the ones present in other quantization procedures. We have discussed this term 
in connection with the Dirac's quantization procedure arguing that its presence is 
connected to the (physical) way of constraining the dynamics on a manifold. 

In this paper we have adopted for po and only the domain with periodic bound- 
ary conditions. Actually there is an infinity of subsets of L^(0,27r) where every op- 
erator we have considered is self-adjoint, i.e. those with iIj{2tt) = e'^'^'^^ip{0) where 



(3 G [0,1[. This issue is clearly exposed in [12, 13| and references therein. One can 
regard the gauge transformation with parameter a in Sec. 3. 2 as a similarity trans- 
formation between these subsets of L^. The (potential) freedom in the choice of the 
domain of definition of the operators is contained in this gauge transformation. 

It would be interesting to elucidate the features of this formalism, in the form 
explicitly including the Lagrange multipliers, in connection with the Faddeev and 
Popov functional technique in quantum field theory |l^ . 



The author would like to thank P. Facchi and S. Pascazio for interesting remarks. 
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